Totally Asymmetric Exclusion Process with Hierarchical Long-Range Connections 
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A non-equilibrium particle transport model, the totally asymmetric exclusion process, is studied 
on a one-dimensional lattice with a hierarchy of fixed long-range connections. This model breaks the 
particle-hole symmetry observed on an ordinary one-dimensional lattice and results in a surprisingly 
simple phase diagram, without a maximum-current phase. Numerical simulations of the model with 
open boundary conditions reveal a number of dynamic features and suggest possible applications. 



I. INTRODUCTION 

Physicists have long hoped to understand non- 
equilibrium phenomena as well as they understand equi- 
librium phenomena [l| . Certain non-equilibrium systems 
exist which reach a steady state, yet they do not obey de- 
tailed balance required for any equilibrium. The steady 
states are defined by the dynamics rather than an en- 
ergy function. Exclusion processes are widely studied 
as models of particle transport and were first introduced 
for the kinetics of bio-polymerization on nucleic acid tem- 
plates [2l . They have since been related to otherphenom- 
ena such as surface growth [Si] , traffic flow jj, , and 
the statistics of DNA sequence alignment [3]. The to- 
tally asymmetric exclusion process (TASEP) represents 
a rare example of an exactly solvable model with a non- 
equilibrium steady state that allows a deep, analytic in- 
sight into the phenomenology of critical behavior beyond 
thermodynamic equilibrium d, [l^. TASEP and re- 
lated models exhibit non-equilibrium phase transitions 
which have no analog in equilibrium systems, such as a 
phase transition in one dimension. 

TASEP describes particles conducting nearest neigh- 
bor jumps along one direction within a line of sites. There 
is no passing, and jumps of particles to an occupied for- 
ward site are excluded, leading to jamming. The model 
has been solved first in a mean-field treatment Q and 
subsequently in full detail d, [lOl. and has since inspired 
a large number of variations [ill. [T3 |. Yet, the phase 
diagram has proven quite robust under those changes. 
Even admitting long-range jumps, which allows particles 
to pass forward by a stochastic long-range jump accord- 
ing to a Levy distribution if the target site is unoccupied, 
leaves the phase diagram surprisingly unchanged [13| . 

The purpose of this paper is to introduce jumps into 
the TASEP in a non-stochastic way by using a network 
with predetermined long-distance jumps and to study it's 
effect on the phases and transitions. Such a quenched 
structure is provided by the recently introduced network 
HNS. HNS has a hierarchical structure, combining a sim- 
ple one-dimensional backbone with a sequence of long- 
range links. Interesting properties for other statistical 



models on HNS have already been described in Ref. [Ij] . 
The process on HNS (HNS-TASEP) might prove to be a 
benign enough extension beyond that in one dimension 
(Id- TASEP) such that analytical insights remain possi- 
ble, even though introducing quenched long-range con- 
nections removes the particle-hole symmetry. Here, we 
show numerically that HNS alters the phase diagram sig- 
nificantly. In fact, only two phases remain, separated by 
a sharp first-order transition. Our results suggest that 
there might indeed be a simple solution for this model 
which would take the analytic treatment of exclusion 
processes beyond one dimension. The observed behav- 
ior suggests that the hierarchical lattice used here also 
provides an efficient switch for a simple, one parameter 
storage-and-release system. 

There are some real systems to which an exclusion 
process with particles passing each other applies. For 
instance, as pointed out in Ref. 01 , some proteins regu- 
late genes by binding to DNA and search for a specific 
target site by continuous dissociation and re-association 
with the DNA. In the dissociation process, they may fully 
dissociate or stay within the range of electrostatic forces 
of the DNA. In another application, HNS-TASEP could 
model traffic d, d, 0| with expressways, which also func- 
tion as quenched shortcuts. The hierarchical but geo- 
metric structure of HNS may even apply to a multi-level 
transport system such as a package delivery service with 
many door-to-door vans, a number of intercity trucks and 
trains, and a few transcontinental fiight routes. 

Our discussion proceeds with a review of the ordinary 
Id- TASEP in the next section, followed in Sec. IIIII by a 
description of the hierarchical lattice geometry employed 
here. In Sees. IIVI and |Vl we present our numerical re- 
sults and analytic approach, respectively. In Sec. IVII we 
discuss the implications of our findings and provide our 
conclusions in Sec. IVIII 
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II. TASEP ON A LINE 

The totally asymmetric exclusion process in one di- 
mensional (Id- TASEP), where particles always move in 
one direction, is a model which has several phases with 
first and second-order transitions 8, 9, 10). It is defined 
on a one-dimensional lattice of length L. Each site of the 
lattice, labeled by i, is either occupied or unoccupied by 
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FIG. 1: Phase diagram of the Id-TASEP with high density 
(HD), low density (LD) and maximum current (MC) phases. 
The injection rate is a and the removal rate is /3. 



a particle, and accordingly has an occupation number 
which is either 1 or 0. Particles on the lattice may only 
move in one direction, which we will say is to the right, 
and they may hop one site to the right only if that site 
is unoccupied. Typically, the system is updated random 
sequentially, that is, particles are selected to move one at 
a time in random order. The average density p is defined 
as N/L, where N is the number of particles in the system. 
The average current J is defined as the average number 
of particles that move through a point on the lattice per 
time step. For periodic boundary conditions, N is fixed, 
while for open boundary conditions, as considered here, 
N is allowed to fluctuate, which induces several phases 
characterized by different properties for p and J. The 
two open boundaries are connected to a large reservoir 
of particles, so that the rate of particles hopping onto the 
lattice at the first site is a, and the rate of particles re- 
moved from the last site is /3, providing two independent 
parameters. Both define a parameter space for which the 
different phases appear, see Fig. [T] Their rates are be- 
tween and 1, as no more than one particle may appear 
or disappear on the boundary in a single trial. 

Fig. [T] shows three phases distinguished by their aver- 
age densities and currents in the steady state, first found 
in Ref. \§\. The low density (LD) phase has pld = a 
and Jlo = a(l — a), the high density (HD) phase has 
Phd = 1-/3 and Jhd — /3(1 — /9), and the maximum 
current (MC) phase has pMC — 1/2 and Jmc = 1/4. 
The HD and LD phases are related because of particle- 
hole symmetry. The transitions between the LD and HD 
phases and the MC phase are continuous. The line be- 
tween the LD and HD phases is called the "shock phase" 
(SP). This transition is not continuous and the system 
reaches a state of coexistence between the two phases. 
On the lattice there is a region with low density and a 
region of high density separated by microscopically small 
"shock" , which diffuses between the ends of the lattice. 

Id-TASEP has been solved exactly with recursive 
equations [1, [loj , as well as a more advanced matrix for- 



mulation [9| . The process is completely described by the 
change in occupations on affected sites for an update at 
a bulk-site i during the time interval \t^t + dt]. Such an 
update merely alters site i and i + 1: 

n{t + dt) = T,(t)T,+i(t) (1) 

all other sites remain unchanged. Special treatment ob- 
tains for each of the two boundary sites. When the up- 
date selects to inject a new particle into the system, that 
particle attempts to occupy site i = 1 with probability 
a, if that site is open. When the last site, i ~ L, is 
selected for an update, it unloads an occupying particle 
with probability /3. 

These master equations, as well as those boundary con- 
ditions, can be averaged over noise, eliminating all fluc- 
tuations, and written as single differential equation de- 
scribing the time evolution of the system. Assuming the 
existence of a steady state, dt (r^) = 0, leads to a system 
of algebraic equations at most quadratic in (r^), which 
can be solved recursively. 

Many variations of TASEP have been studied, such as 
those with parallel updates, multiple species of particles, 
and extended particle sizes [l3|. Remarkably, the phase 
diagram is essentially the same for many variations. For 
example, one variation studied recently introduced long 
range jumps so that particles may jump past each other 
(if the target site is unoccupied) to study the changes to 
the phases and transitions [l3|. These jumps are stochas- 
tic Levy flights, i. e. a distance I is reached with probabil- 
ity pi ^ depending on a parameter a > 0. Even 
with long jumps and reordering ("passing") of particles, 
the diagram has the same three phases as in Fig. [T] 

III. TASEP ON THE HIERARCHICAL 
NETWORK HNS 

HNS is a network with a fractional dimension, which 
was introduced with the intention of studying small- 
world phenomena analytically [14] . Small- world phenom- 
ena are found in many natural and man-made systems, 
such as neural networks and the internet [lB|. The net- 
works are characterized by having a mixed structure , 
with connections between geometrically nearby neighbors 
as well as long distance connections, which drastically re- 
duce the typical path length between any two nodes. The 
HNS network does not have the same mean-field proper- 
ties usually associated with such networks. For instance, 
average path lengths scale with VL instead of ln(L) with 
system size L. But HNS is constructed in a hierarchical 
manner that is conducive for the renormalization group, 
which is a technique that takes advantage of the self- 
similarity of systems [l^ . 

HNS consists of a one dimensional line as a backbone 
with L = 2'^ + 1 sites, where fe is a positive integer which 
defines the number of hierarchies. To make the long dis- 
tance connections, we consider an integer i < k, which 
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FIG. 2: The HNS network consists of a one-dimensional back- 
bone and a hierarchy of long distance connections. 



defines the level of the hierarchy, and j, an integer which 
parameterizes the connection in a hierarchy. All sites 
(except for n = 0) are then uniquely represented by 



n = 2\2j + l). 



(2) 



For example, for i = 0, n runs over all the odd integers, 
and 4 = 1 makes n to be all integers once divisible by 2 
(i.e. 2, 6, 10, ...), etc. Connections are made between 
neighbors within the hierarchy, so for z = 0, site 1 con- 
nects to 3, 5 to 7, etc, and for i = 1, 2 to 6, 10 to 14, etc. 
It is possible to make a more complex network, HN4, by 
connecting two neighbors in the hierarchy to every site. 
Fig. m shows the HNS network for A; = 5. The first, mid- 
dle, and last site are special in that they have no long 
range connections (the first and last site, and L, are 
not shown). All the other sites are connected to three 
other sites. The distance between the ends of the HN3 
network are known to be proportional to a/Z, which is 
similar to the diagonal on a square lattice with L sites. 

In an implementation of TASEP on the HNS network 
(which we shall refer to as IIN3- TASEP), on sites with 
a long-range forward link, particles have the possibility 
to move to two different sites. To obtain an interesting 
dynamics we decided to have the chosen particle attempt 
a jump to the long distance site first, and if that site is 
occupied, to reach the short distance site as in Id- TASEP. 
For such a site, representing half of all sites, jamming is 
somewhat alleviated, as it can free itself with a higher 
probability. In turn, the other half of all sites, possessing 
one extra incoming link instead, occupation and jamming 
is far more likely. 

We will also study a family of models that interpolates 
between Id-TASEP and HN3-TASEP, using a probabilis- 
tic update rule: At each update, with a probability r, a 
long-range jump is attempted first but if such a jump 
does not succeed, a nearest-neighbor forward-jump is at- 
tempted. For r = 0, no long-range jump ever occurs, rep- 
resenting the Id- TASEP case, while r = 1 corresponds 
to the HN3-TASEP case. We find that the nature of 
TASEP changes discontinuously for r > 0, numerically 
signaled by the disappearance of the MC phase (and the 
associated 2n(i-order transitions), while the shock per- 
sists at any r. We suspect that such discontinuity can be 
attributed to the broken particle-hole symmetry, which 
only holds for strictly r = 0. 





FIG. 3: to compare with that of Id-TASEP in Fig. [T] The 
density p (top) and current J (bottom) of TASEP on HNS on 
a grid of 25 values each for a and (3. 



IV. NUMERICAL RESULTS 

We implemented a Monte Carlo simulation for HN3- 
TASEP as follows. First a particle is randomly chosen 
from N particles on the lattice plus one virtual parti- 
cle. If a lattice particle is selected, it is moved forward 
one site if that site is unoccupied. If a virtual particle is 
selected, a random number between zero and one is gen- 
erated, and if it is less than a, a particle is placed at the 
first site of the lattice unless it is occupied. If a particle is 
selected which is at the end of the lattice, a random num- 
ber is generated, and if it is less than P, that particle is 
removed. A sequence of A^-|- 1 attempts to move particles 
constitute one Monte Carlo Sweep (MCS). The number 
of lattice sites for every result in this report was set at 
L — 1023, unless otherwise stated. The simulations were 
run for 10^ MCS, and the first 10^ MCS were discarded 
to allow the system reached a steady state. Current and 
density for the system were recorded every 100 MCS and 
averaged over to characterize the steady state. 

The behavior of HN3-TASEP is quite different from Id- 
TASEP [8] , which can be seen when p and J are plotted 
for a grid of 25 by 25 values of a and /3 in Fig. [3] The 
HD and LD phases are present, but the MC phase is 
conspicuously missing. The magnitude of p and J are 
significantly altered. (Note that the current is measured 
here at the exit site, which every particle moving through 
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FIG. 4: Phase diagram of HN3-TASEP with only a high den- 
sity (HD) and low density (LD) phases, to compare with that 
of Id-TASEP in Fig. [1] The data points are extracted from 
Fig- El the dashed line represents a simple polynomial fit, and 
the dotted line is crude approximation to the phase line re- 
sulting from the linearized mean-field equations in Sec. [V] 



the system must pass; the same is not true for most other 
sites!) The density is much higher in the HD phase and 
much lower in the LD phase when the HNS shortcuts 
are present. In fact, throughout each phase, the density 
remains almost constant. The lattice is able to fill itself 
with particles more efficiently in HD and remove them 
more efficiently in LD when there are more connections 
between sites. Based on the density plot in Fig. [31 in 
Fig. m we extract the outlines of a phase diagram for HN3- 
TASEP to compare with that of Id-TASEP in Fig. [2 

Considering that HNS is a hierarchical network, it is 
worth noting that the results are smooth and not hetero- 
geneous in any complicated fashion (unlike the average 
occupation on sites, see below). Similar to Id-TASEP, in 
the HD phase the current does not appear to vary with 
a, and in the LD phase it does not change appreciably 
with (3. To make sure the transition is first-order, p is 
plotted for values of a = 1 and /3 between 0.7 and 0.8 
for two lattice sizes in Figs. |5l The transition is less pro- 
nounced on a smaller lattice, which indicates the transi- 
tion is most likely sharp in the thermodynamic limit: The 
data indicates that in the thermodynamic limit L ^ oo 
the lattice is completely filled almost everywhere in the 
high-density phase, and it is completely empty in the low- 
density phase. This result enables a simplifying Ansatz 
to the equations in Sec. FVl below. 

If one considers a single particle moving through an 
empty Id lattice, it moves one site every step and u = 1. 
Since in HNS the shortest end-to-end path is of length 
~ Vl [i3| , the time for a particle to traverse it is ~ -s/L 
and, hence, v = L/T w \/L. The velocity is dependent 
on the lattice size, and is unbounded for large lattices, 
although this definition of velocity is based on the notion 
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FIG. 5: Top: Density during the transition from HD to LD for 
lattices of size L = 255, 511, 1023, and 2047 at q = 1. The 
total number of Monte Carlo Sweeps (MCS) for each data 
point was 5 lO". Bottom: The density for the different values 
of f3 plotted against 1/L. The values range over f5 = 0.7 . . . 0.8 
in steps of 0.08. The dashed line in the middle corresponds 
to /3 = 0.748, which appears closest to the transition, exhibits 
a significant error as the stationary "shock" state is hard to 
sample. But almost all data below extrapolates to a value of 
p consistent with full packing, while all data for values of /3 
above the transition extrapolates to an empty lattice, p ~ 0. 



that the length of a lattice is equal to it's number of sites. 
In the LD phase the particles can follow the shortcuts, 
since they are usually unoccupied. In the HD phase most 
sites, including ones at the end of shortcuts, are occupied, 
so the particles have few opportunities to take a shortcut 
and they travel along the backbone, which limits the ve- 
locity to V < 1. Fig. |S1 shows the velocity for a parameter 
range across the transition. In the HD phase the veloc- 
ity is indeed < 1, but in the LD phase it is ^ 1. The 
velocity can also be examined for the case (3 = 1, where 
the system is always in the LD phase, as shown in Fig. [71 
For small a, the number of particles is very low and they 
can take all of the shortcuts. The velocity is on the order 
of 16 « yi/2. As a increases, some of the shortcuts 
are occupied and the particles take trajectories that are 
a mix between shortcut and backbone movements, which 
decreases the average velocity. The -v/X-dependence of 



5 




0.2 0.4 0.6 0.8 1 



P 

FIG. 6: Average velocity of particles for a = 1 and varying /3. 
The velocity is n < 1 in the HD phase {(3 < 0.75), and n > 1 
in the LD phase. 
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FIG. 7: Average velocity of TASEP on HNS for /? = 1. 



the transit time can be demonstrated through finite size 
scaling. Fig.[5]shows the average transit time at a low in- 
jection rate and high removal rate for many lattice sizes, 
and it approximately follows ^/L. The stair-casing effect 
between even and odd values of the hierarchy index k is 
due to the fact that HNS is strictly self-similar only for 
every second recursion of the hierarchy. In particular, 
as particles move left to right, they take every shortcut 
available, and on networks with odd k they miss the op- 
portunity to take the largest shortcut, which is part of 
the shortest path available. 

One may also keep track of the number of time steps it 
takes for each particle to cross the lattice. Fig. [5] shows 
a histogram of the transit times for f3 —1 and a small 
value of a, so the particles are removed at the maximum 
rate and injected at a slow rate. There is a bimodal dis- 
tribution which suggests particles are only taking a few 
distinct paths through the network, such as the shortest 
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FIG. 8: Average transit time versus lattice size L for a — 0.01 
and 13—1. The transit time shows a -^/L dependence, which 
means particles are taking the shortest path on HNS. The 
stair-casing is due to particles on networks with odd k only 
taking the second shortest path. 
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FIG. 9: Distribution of the transit times at a = 0.1, and 

/3= 1. 



and second shortest paths. It is interesting that the bi- 
modal distribution does not show up for other parts of 
the parameter space. When a is smaller, the particles 
follow only the shortest path, and there is only one peak. 
The width of the peak originates with the update pro- 
cedure, as not every particle is updated at every Monte 
Carlo sweep, thus leading to fluctuations in transit times. 
When a becomes larger there is some blockage of short- 
cuts and the particles immediately follow many different 
paths. There is only one peak in the distribution of tran- 
sit times with a long tail and there is no clear separation 
between fast moving particles and slow moving particles. 
The number of particles taking the shortest path is prob- 
ably very small for any set of parameters except when the 
density is very low. 
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a.= 1 , |3=D.8 averaged over 1DDDD steps 




FIG. 10: Time evolution of the shock phase in TASEP on HNS 
on a lattice of size L — 255. Each white pixel is an occupied 
site. Time is represented by 100 Monte Carlo sweeps. 



We note that the shock phase looks very different on 
HNS from that in Id- TASEP. Fig. [TUl shows a typical time 
evolution of the occupation on the lattice at a point in 
the parameter space where p was roughly 0.5 and can 



be compared to the shock phase in Id- TASEP [ll[. It is 
clear that there is coexistence of the LD and HD phases, 
but there appears to be a hierarchy of boundary sites, 
depending on the range of the long connections of the 
originating site, as we might expect that there is more 
movement near sites with the longer connection. 

A similar picture about the intricate internal dynamics 
within the lattice emerges from the average steady-state 
occupation {n) for each site, shown in Fig.[TTJ Smoothed 
over many sites, one could argue that the general trend 
is similar to that in Id- TASEP, where in LD the bulk 
density is constant, cumulating in a defined boundary 
layer near the exit, while in HD there is such a layer 
at the entrance, followed by a constant density plateau 
in the bulk [1]. Although in HN3-TASEP correspond- 
ing boundary layers are visible, the site-to-site density 
is extremely rough, owing to the heterogeneous mix of 
incoming and outgoing long-range jumps that belong to 
very different levels of the hierarchy and thus have very 
different efficiencies in transmitting particles. Our linear 
approximation to the steady-state equations in Sec. |V| 
manages to reproduce these heterogeneities very well. 



V. ANALYTIC TREATMENT 

To obtain a set of master equations such as those in 
Eq. H]) for HN3-TASEP, we have to distinguish two dif- 
ferent kinds of sites. An IN-site i simply has a long-range 
connection to a site i — I, preceding it by a distance I in 
the lattice, in addition to its immediate predecessor and 
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FIG. 11: Average occupation (r^) in the steady state on sites 
along the lattice backbone in HNS- TASEP in LD (top) and 
HD (bottom). In low density (LD), sites with the longer- 
range connections show more occupation as way-stations of 
short paths through the lattice that bypass particles away 
from sites belonging to lower levels of the hierarchy. In turn, 
in high density (HD) low-level sites get jammed up while sites 
with longer-range connections (especially forward) have an 
easier time to empty themselves. Note the boundary layers of 
increased jamming near the exit in LD and that of depletion 
at the entrance in HD. 



successor sites i — l and i + l. Generically, at any level in 
the hierarchy, updates at an IN-site i do not affect either 
of the predecessor sites, hence, the equations are similar 
to those in Eqs. ([T]): 

nit + dt) = n{t)n+i(t), (3) 
Ti+i{t + dt) = n+i{t) + [i^n+i{t)]n{t). 

In turn, when the update occurs at an OUT-site i with 
a long-range link to a successor site i + l, that site 
the immediate successor site i -I- 1, and the long-range 
successor site i -\- I are affected in a novel way in HN3- 
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FIG. 12: Average steady-state occupation (r^) ~ 1 — on 
HN3-TASEP, as obtained analytically in the linearized ap- 
proximation (ti <^ 1) in Eq. p5|) for HD. Aside from an arbi- 
trary overall scale for the e^, the pattern corresponds in great 
detail to that in Fig. 1111 except for the boundaries. 



TASEP. In our choice of moving particles preferentially 
along the long links leads to a three-stage update process: 
Only if the long jump is blocked, a short jump to the 
successor site is attempted; if that is blocked as well, the 
particles remains on its site. These choices are expressed 
through the following equations for the updated site, its 
successor, and its long-range forward neighbor: 

n{t + dt) = nit)T,+i{t)T,+i{t), (4) 
n+i{t + dt) = Ti+i{t) + Ti(t)[i-n+i{t)]Ti+i{t), 

Ti+l{t + dt) = Ti+l{t)+Ti{t)[l-n+l{t)]. 

Note that these equations are inherently cubic in the 
dynamic variables. Another complication is the lack of 
translational invariance, as these equations depend on a 
forward-distance which itself depends strongly on the 
hierarchical level that site i belongs to. At least, the 
boundary conditions, affecting sites i — Q and i — L, are 
identical to Id-TASEP. 

To make any progress at all, we already at this point 
consider the mean- field limit. The mean- field limit av- 
erages the dynamic variables over the noise, eliminat- 
ing fluctuations and correlations, i. e. we set (tiTj) ~ 
(Ti) {tj) = TiTj, and allows us to arrive at a set of rate 
equations for the continuous variables Ti . We have to take 
full account of the lack of translational invariance along 
the line in HNS by addressing the hierarchical level any 
site belongs to. All sites on the lowest level are of odd 
index 2j + 1, and they are alternately IN and OUT-sites. 
Say, all sites 4j -|- 1 are OUT-sites. While it is then clear 
that it has two successors, one and two steps ahead, and 
a single predecessor site, the latter itself may be an IN 
or an OUT-site. But depending on that, its update will 
affect our Aj + 1 site differently, and we make the further 
simplifying assumption of an equal balance between both 



possibilities. Hence, on average, site 4j -I- 1 is changed as 

Ti^ + l{t + dt) = {1- dt)Ti^ + l{t) (5) 
+dt Tij+l (t)Tij+2 {t)Ti]+3 (t) 

+ y^4j(<) [l-T4, + l(t)] 

where the distance I furthermore depends on that pre- 
decessor site at 4j. In order, the terms in Eq. ([5]) refer 
to either nothing changing with probability 1 — dt, our 
reference site 4j -I- 1 being updated itself [see the first 
of Eqs. (H])] with probability dt, or the predecessor site 
4j being update with probability dt/2 for each of the 
two scenarios given. Similar considerations holds for the 
IN-sites at 4j + 3, with yet another complication aris- 
ing from the incoming long-range predecessor site 4j + 1, 
whose potential update adds another term: 

Tij+3it + dt) = il-dt)T4j+3{t) (6) 
+ dtT4j+3it)T4j+i{t) 

+ yT4j+2(t) [l-r4j+3(t)] 

+ y Tij+2it)Tij+2+l{t) [1 - Tij+sit)] 
+dtT4j + i{t) [1 - T4j+3(t)] . 

At all other levels in the hierarchy, matters somewhat 
simplify, as any even site clearly has an odd-indexed site 
preceding and following it, making the effect of long- 
range bonds fully apparent. In particular, sites 2{2j + 1) 
at the next-to-bottom level all have an odd OUT-site 
4j -I- 1 preceding it and the corresponding IN-site 4j -I- 3 
following it, independent of whether they themselves are 
OUT or IN-sites. Still, in their own update behavior, 
OUT-sites at 2(4j + 1) and IN-sites at 2(4j + 3) at this 
level differ, and we get 

T2(4j + 1) {t + dt) = (1 - dt) r2(4j + i) (t) (7) 
-l-dt T2(4j + i) (t)r2(4j + i) + i(t)r2(4j+3) (i) 
-Kdi T2(4j + i)_i(<)T2(4j+;) + i(i) [l - T2(^4j+l){t)] , 

and 

T-2(4j+3) it + dt) = {l- dt) T2(4j+3) {t) (8) 
-l-di T2(4j+3)(t)T2(4j+3) + i(i) 

+dtT2(4:j+3)-l{t)T2{ij+3) + l{t) [l - T2(4j+3)(0] 
+dtT2(4j + l)(i) [1 - T2(4j+3) W] ■ 

Finally, at any higher level « > 2 of the hierarchy, all 
sites 2* (2j -I- 1) only have odd-indexed IN-sites as prede- 
cessor and odd-indexed OUT-sites as successor. We get 
for OUT and IN-sites, respectively: 

T2i{4j + l){t + dt) = (1 - dt) r2i(4j + i)(t) (9) 
+ T2.(4j + i)(t)T2»(4j + i) + l(t)T2.(4j+3)(t) 
-fdtr2.(4j + i)_i(i) [1 - T2.(4j + i)(i)] , 
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and 



Special rules apply for the boundary sites in Eqs. (fT4|) : 



T-2.(4j+3) {t + dt) = (1 - dt) T2.(4j+3) (t) (10) 
+dt T 2» (4j+3) 2' (4j+3) + l (i) 
+diT2.(4j+3)_i(t) [l - T2.(4j+3)(i)] 
+diT2.(4j + l)(i) [1 - T2.(4j+3) W] • 

In the steady state, we consider t ~ i + rft — > oo. Then, 
the above equations somewhat simplify, and we get from 
Eqs. dSlini): 

Tn = T„r„+ir„+2 + -r„_i (1 - r„) (1 + T„_i+;) , (11) 

7"™ = TmTm+i + -rm_i (1 — Tm) (1 + rm_i+i) 
+'rm-2 (1 - Tm) , 

at the lowest level, where we have set n — 4j + 1 for 
OUT-sites and m = 4j + 3 for IN-sites. Similarly, from 
Eqs. dZKHl), we get 

Tn = TnT„ + lTn+4 + f/i-l (1 - T„) T„+i, (12) 
Tm — Tm^m+l ^" '^m — 1 (1 '^m) 
+r„i_4 (1 - T,„) , 

here setting n = 2(4j + 1) for the OUT-sites and m — 
2(4j-|-3) for the IN-sites. Finally, the general case in the 
hierarchy from Eqs. ([ WTU)) reduces to 



Tn = TnTn+lTn+2^ + 1 + Tn-1 {I — Tn) ■ 
Tm — TniTni^l ^ Tni—l (1 '7'm) 
+7'm-2^+i (1 - Tm) , 



(13) 



with n = 2'(4j + l) for the OUT-sites and m = 2*(4j + 3) 
for the IN-sites for i > 2. Since the first and last bond in 
the lattice do not get bridged by a long-range jump, see 
Fig. [21 the boundary conditions are similar to those for 
Id-TASEP fSl: 



To (1 - Ti) = a{l- To) , 

f^TL = TL-1 (1 - Tl) , 



(14) 



with L = 2'=. 

Although the steady state has simplified the equations 
somewhat, their main difficulties remain their lack of 
symmetry and their cubic order. Yet, our numerical re- 
sults in SeclIV) particularly Fig.[3l suggest a very simple 
phase structure with ^ 1 — e(a, (3) in the bulk through- 
out HD, and 6 {a, (3) throughout LD, with e,(S <C 1. 
Already a linear approximation in HD provides some illu- 
minating insight. |18j To first order in e^, the distinctions 
between the different orders in the hierarchy, as expressed 
in Eqs. (|imi3p . disappear to leave for all < z < fc — 2 
and < j < 2'=-2-*: 



e2'(4j + l) ^ e2'(4j + l) + l + e2'(4j+3)) 

1 

e2'(4j-|-3) ^ 2'^2*(4j+3) + l- 



(15) 



(l-eo)ei 
/3(l-eL) 



aeo, 

(1 - CL-l) CL, 



(16) 



which we leave unexpanded for now, and the central site: 



e2fc-i 



£2^-1-1, 



(17) 



each of which lacks a long-range bond. The solution 
for all Ci, while not presentable in closed form, is eas- 
ily obtained in 0{L) steps in terms of e^, starting from 
the exit backwards. Fig. [T^] demonstrates the quality of 
this approximation in accounting in great detail for the 
Weierstrass-like hierarchical nature [l3| of the site densi- 
ties in HD found in Fig. 1111 In particular, the pattern on 
bulk sites is well represented, aside from an overall scale 
represented by the value of that depends on a and (3. 
Naturally, sites near the boundaries are less- well approx- 
imated, since the are not sufficiently small there. 
Current conservation on the first and last bond dictates 



£1 



(18) 



which is automatically satisfied by Eqs. (jlSp . Together 
with Eq. (flB)) . Eq. ([18]) relates the boundary conditions 
together, which allows to fix the arbitrary scale in terms 
of a and (3. More interesting, we can estimate the critical 
line {ac, /3c) between HD and LD as the location where 
the linear approximation breaks down. Eliminating the 
other unknowns, we find the nontrivial relation involving 
both boundary conditions: 



aeo 



1 



1 aeo 

2 1-eo 



1-eo 1 



1-eo 



(19) 



that fixes eo{a,f3). In turn, any such relation is only 
reasonable for sufficiently small eg. Since eq varies be- 
tween and 1, either extreme providing a trivial or singu- 
lar result, respectively, we chose (somewhat arbitrarily) 
a generic in-between value of eo = 1/2 to estimate the 
phase line. The resulting relation, a — 1 + P — y^l + (3-^, 
is also indicated in the phase diagram in Fig. [Jj 



VI. DISCUSSION 

We found that there are many striking differences be- 
tween Id-TASEP and TASEP on HN3. It is obvious that 
particle-hole symmetry is broken from the asymmetric 
phase diagram, where the phase boundary between HD 
and LD curves and favors a larger LD phase. The sym- 
metry is broken due to the particle's preference to take 
shortcuts, and Fig. [THl illustrates the most obvious case 
of this symmetry breaking. A particle may only move 
from site 1 to site 3, but a hole may not move from site 
2 to site 1 since the particle may not move from site 1 to 
site 2. A full lattice with a single hole moving to the left 
cannot take every shortcut as does a single particle on an 
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1 2 3 



FIG. 13: Particle-hole symmetry is broken by the particle's 
preference to take the shortcut. A particle may move from 
site 1 to site 3, but a hole may not move from site 2 to site 1. 



0.5 




qI ^ ^ ^ ^ 1 
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FIG. 14: Plot of the current J(a) for fixed /3 = 1 in the one- 
parameter family of models described in the text for < r < 1 
in a system of size L — 1024. For r — only, the case of 
Id-TASEP, particle-hole symmetry is obeyed and phase tran- 
sition into an MC phase occurs. For all r > 0, the symmetry 
is broken and J{a) varies smoothly, similar to r = 1, the case 
of HN3-TASEP. 



empty lattice moving to the right. The model favors the 
efficient flow of particles, and this is why the LD phase 
is larger and free flow persists for some values of a > /?. 

The importance of this symmetry is highlighted by de- 
signing a one-parameter family of models that interprets 
between Id-TASEP and HN3-TASEP. In these models, 
particles attempt a long-range jump with probability r, 
< r < 1. Should no jump result, either because no at- 
tempt was made or because the attempt failed, a nearest- 
neighbor jump is attempted. For r = 0, only nearest- 
neighbor jumps occur, corresponding to Id-TASEP. In 
turn, for r = 1, a long-range jump is always attempted 
first, as in HN3-TASEP. We can characterize the behav- 
ior of this model with r sufficiently by presenting our 
numerical results for the current J as a function of a for 
fixed (3 = 1. Fig. [T3] shows how the smooth variation 




FIG. 15: Density p (top) and current J (bottom) of TASEP 
on HN3 with particles having 50% chance of only attempting 
long distance jump and 50% chance of only attempting jumps 
along the backbone. 



of J{a) found for HN3-TASEP in Fig. [3] evolves from 
r = 1 towards the 2nd-order phase transition between 
LD and MC phases at r = 0. Although finite-size effects 
(at L = 1024 here) obscure the bend in J at a = 1/2, it 
is quite clear that only for strictly r there is a MC 
phase with a constant-current plateau. For any r > 0, 
J(a) appears to remain a smooth function without the 
emergence of a plateau. Of course, for all such r > 0, 
particle-hole symmetry is broken. 

We can formulate an update rule on HNS that does 
preserve particle-hole symmetry. For instance, the parti- 
cle chosen for an update could attempt a long-range or 
a nearest-neighbor jump with probability r and 1 — r, 
respectively, but not explore the alternative if such an 
attempt fails. Unfortunately, although r = again corre- 
sponds to Id-TASEP, r = 1 does not attain HN3-TASEP 
but another version of Id- TASEP restricted to ^ VL 
sites only. In Fig. [15] we display the density and the cur- 
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rent as a function of a and /? for the case r = 1/2. The 
symmetry with respect to a = /? is easily visible. As 
in Id-TASEP, there are HD, LD, and MC phases with a 
shock hne between HD and LD. The most notable differ- 
ences are that J > 1/4 in MC and the steep transitions 
in the density at each phase boundary. 

Returning to HN3-TASEP, we found in the LD phase 
that the particles move freely across the lattice. As the 
injection rate a is increased, the density and current in- 
crease, and there are few collisions, but the network is 
still able to transport particles quickly. At the transition 
into the HD phase, jams snowball throughout the system 
and soon fill in the whole lattice. The shortcuts are likely 
to be blocked off and the particles slowly crawl along the 
backbone. One may make an analogy to vehicular traf- 
fic. Cars like particles can move along a one lane road, 
and they may not pass each other. They appear at the 
beginning of the road at a rate a, and disappear at the 
end of the road at a rate /?. At some density of cars the 
average velocity drops significantly and the whole road 
is jammed. If we add highways then cars may enter a 
highway, go faster, and when they get off they end up 
in front of their peers. In free fiowing traffic, highways 
reduce the amount of time for cars to reach their desti- 
nation and lower the density, given the same boundary 
conditions. Traffic jams still occur with the addition of 
highways, and there are more possibilities for cars to ma- 
neuver around each other and fill in gaps so the density 
can be very high. 

An interesting property of TASEP on HNS is that the 
network uses mainly it's long range connections to trans- 
port particles, and its backbone to 'store' them. In the 
LD phase the particles flow freely through the highway. 
As a is increased, particles may be pushed into the less 
visited backbone sites, and the current and density in- 
crease. An inverse process happens in the HD phase, 
where as is increased more particles can be pulled from 
storage, the current goes up, and the density goes down. 
There is no possibility for the current and density to stay 



constant as both parameters are changed simultaneously, 
hence, there is no MC phase on HNS. 

This storage mechanism could conceivably be useful 
in a real system. The transition between the HD and 
LD phases changes the system from very high density 
to very low density, while the current does not change 
significantly. Only a small local change in the rate of 
particle injection or removal at the boundaries can induce 
a dramatic global change in the number of particles. We 
have measure the effect of a protocol whereby at a = 1 
the value of (3 is switched between 0.7 (HD phase) and 
0.8 (LD phase), see also Fig. [5] Unfortunately, it takes a 
long time, about a factor of 1000 longer than the typical 
transit time of a particle, to squeeze the excess particles 
through the single exit site, empty out the lattice and re- 
establish the steady state at /? = 0.8. In turn, jamming 
up the system by re-setting to (i — 0.7 attains it steady 
state at least an order of magnitude faster. 

VII. CONCLUSIONS 

In this work we found that a variation of the TASEP 
with fixed long distance jumps can lead to significant 
changes in its phases, notably, the disappearance of the 
maximum current phase. This is surprising since the 
TASEP phases are considered to be robust to many 
changes. These changes were described qualitatively as 
a result of the particles velocity becoming essentially un- 
bounded along paths utilizing the long jumps. 

Despite the increased complexity of the rate equations 
for this model, the much-simplified phase structure we 
found allowed the possibility of an analytic approach 
that, even if approximate, should provide a good qualita- 
tive and quantitative description. While a more detailed 
solution eludes us here, it remains a worthwhile goal for 
the future, as it would provide novel insight into and 
control over a process that has stimulated significant ad- 
vances in the understanding of non-equilibrium critical 
phenomena. 
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